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Abstract. Maximum electric fields of Langmuir waves 
at planetary foreshocks are estimated from the thresh- 
old for electrostatic decay, assuming it saturates beam- 
driven growth, and incorporating heliospheric variation 
of plasma density and temperature. Comparison with 
spacecraft observations yields good quantitative agree- 
ment. Observations in type III radio sources are also 
in accord with this interpretation. A single mechanism 
can thus account for the highest fields of beam driven 
waves in both contexts. 

1. INTRODUCTION 

Electrostatic oscillations near the plasma frequency 
have been observed in the foreshocks of Venus, Earth, 
Mars, Jupiter, Saturn, Uranus, and Neptune [ Gur - 
nett et al., 1989; Grard et al., 1991; Cairns and Gar- 
nett, 1992; Hospodarsky et al., 1994; and the refer- 
ences cited therein]. These oscillations have generally 
been interpreted as Langmuir waves. The observed field 
strengths systematically decrease with heliospheric dis- 
tance [Cairns and Gurnett, 1992], as shown in Fig. 1. 

The currently accepted model for production of Lang- 
muir waves upstream of planetary bowshocks is that 
they are generated via instabilities of electron beams 
that are ejected from the bowshock and propagate along 
the interplanetary magnetic field into the foreshock. Be- 
cause of this similarity and the fact that all foreshock 
plasmas are weakly magnetized (plasma frequency u p 
electron cyclotron frequency) we are encouraged to 
seek an explanation of the trend seen in Fig. 1 in terms 
of a single underlying saturation mechanism, which lim- 
its the fields to below some threshold. 

Foreshock plasmas are very similar to those in the 
solar wind, the main difference being the presence of 
beams generated at the nearby bowshock. Recent work 
on type III solar radio bursts, also generated by beams, 
has provided strong evidence that growth of the most 
intense Langmuir waves L is saturated by the electro- 
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static decay instability L -A- V + S in which a Langmuir 
wave decays into another Langmuir wave V and an ion 
sound wave S [Robinson et al., 1993a, b; Cairns and 
Robinson, 1992a, 1994]; growth of less intense waves 
is retarded by stochastic effects. The product wave L‘ 
is then damped, removing energy from the Langmuir 
population. Evidence exists that this mechanism is im- 
portant in foreshocks [Cairns and Melrose , 1985; Cairns 
and Robinson, 1992a; Thiessen and Kellogg, 1993]. The 
purpose of this paper is to predict the maximum Lang- 
muir fields in planetary foreshocks subject to saturation 
by this mechanism and to compare them with the obser- 
vations shown in Fig. 1. In Sec. 2 we briefly recapitulate 
the theory. A brief summary of the observations is given 
in Sec. 3, and it is shown that their upper bounds agree 
well with theory. 

2. THEORY 

The electrostatic decay process L — > V + 5 has been 
studied by a number of authors, including Robinson et 
al. [1993a], who found the following threshold for the 
process to have a positive growth rate: 

- 3m e V 4 u p 7 z /A^ 

L0 27re 2 7 vs v b v b ’ 

3 [ m 7 L' Av b V ATL ^ 

— 7; Tjo \ Nksle- (J) 

27re 0 7 , V 2 y m e p v b v b 

These equations give the threshold L-wave electric field 
Elo in terms of the electron and ion masses m e and m,, 
the electron and ion temperatures T e and Ti, the elec- 
tron thermal speed V = {ksTe/nie) 1 / 2 , the ion sound 
speed vs = (yksTe/mi ) 1 ^ 2 , 7 = 1 + 3Ti/T e , the beam 
speed v b and its spread Av b , the damping rate 7 l> of 
the waves L' , and the electron number density N . They 
also assume that the product sound waves are relatively 
strongly damped. It should be noted that, although 
Ei 0 is the threshold for electrostatic decay to proceed, 
Langmuir fields can be driven slightly past this level be- 
fore nonlinear damping is high enough to saturate lin- 
ear growth. Robinson et al. [1993a] argued that fields of 
order SElo could momentarily be attained in this way. 
Nonetheless, Elo provides an estimate of the character- 
istic maximum fields attainable, subject to electrostatic 
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decay. Note that Robinson et al.’s [1993a] dynamical 
equations were derived assuming only a single decay, 
but the threshold (1) is independent of this assumption. 

The variation of El q with heliospheric distance r is 
the main focus of this paper. By making some assump- 
tions about the factors in (2), we can estimate the form 
of this variation. Specifically, our assumptions are (i) 
Vb oc V for the beams responsible for the strongest fields 
(i.e., not necessarily all beams), (ii) Avb/vb is constant 
for these beams, and (iii) 7 l< oc u p . These assumptions 
are all plausible if, indeed, there is a single mechanism 
behind the generation of foreshock Langmuir waves. 
Equivalently, one could simply assume that the elec- 
tron velocity distribution has a fixed functional form, 
scaled only by the thermal speed. A third equivalent 
assumption, would be that the electric energy density 
at threshold is a fixed fraction of the thermal energy 
density (cf., Gurnett et al. [1980]). 

The variation of N and T e in the solar wind has been 
studied by a number of authors. At large r, one has 

N oc r~ 2 , (3) 

T c <x r-«, (4) 


with significant ranges in £ predicted theoretically (2/7 
- 4/3) and obtained experimentally (0.24 — 0.70). Scud- 
der and Olbert [1979] and Sittler and Scudder [1980] 
found that a theoretical value of £ = 1/3 fitted the ob- 
servations better than other alternatives between 0.45 
and 4.76 AU; their empirical fit yielded £ = 0.35 ±0.06. 
Phillips et al. [1993] find £ ~ 0.70 from 1.2 to 3.8 AU 
and £ ~ 0.24 from 3.8 to 5.3 AU. Using the above as- 
sumptions, Eqs (2) to (4) yield 


Elo ~ 2.4 


( lL' ( r ^ — 1 — C/ 2 

Vl s — 1 ) VI AU/ 


mV m 1 , (5) 


their frequency does not coincide with the center of one 
of the instrumental channels, the signals are reduced 
by filter roll-off effects [Robinson et al ., 1993b]. Sec- 
ond, the logarithmic compressors used in the ISEE and 
Voyager [W. S. Kurth, personal communication, 1994] 
instruments have exponential time constants ~ 50 ms, 
so shorter signals are underestimated. Third, the high- 
est fields in the various foreshocks may not be among 
those actually sampled by the Voyager spacecraft in the 
short periods of observation (typically 1 hour) close to 
the tangent field line and upstream from the planet (i.e., 
before the beam loses considerable energy). 

The sloping line in Fig. 1, given by (5) for £ = 0.35 
(see below for further discussion) and 7 l> = I0~ 3 u p , 
is consistent with the observed upper bounds. This 
supports our invocation of electrostatic decay as the 
saturation mechanism. The selection effects noted in 
the previous paragraph imply that the observed up- 
per bounds underestimate the true ones, but the de- 
gree of underestimation is likely to be similar in all the 
cases represented by vertical bars, because very similar 
instruments were used throughout. Hence, the slope 
of the heliospheric variation in Fig. 1 should not be 
strongly affected by selection effects. Another source 
of uncertainty is the use of typical parameters at 1 AU 
to evaluate the numerical coefficient in (5) ; variation of 
this normalization by a factor of 3 - 5 is quite possible, 
depending on the solar wind density and temperatures 
and the value of 7 l‘ • 

Phobos 2 observations in the Martian foreshock [Gra- 
rd et al, 1991] yield maximum Langmuir fields of 10 mV 
m _1 , in good agreement with the sloping line in Fig. 1 
(where upper bounds obtained other than by ISEE and 
Voyager are shown as dots). Maximum fields of ~ 1 
mV m -1 were observed in the Venusian foreshock by 


where we have normalized to typical parameters in 
the Earth’s foreshock T e = 1.5 x 10 5 K, 7 = 1.8, 
u p = 1.5 x 10 5 s -1 , Vb = 5V, and A Vb/vb = 0.3. We 
take 7 l> ^ 10 -3 o> p nominally, higher than the values 
of ~ 10 ~ 4 u p seen in the solar wind because the beams 
are slower here, leading to more damping by the back- 
ground thermal plasma. Conversely, 7 l> is lower than 
the ~ 10 ~ 2 u p found for beams deeper in the Earth’s 
foreshock, which are slower still [ Newman , 1985; Fitzen- 
reiter et al., 1984]. 

3. COMPARISON WITH 
OBSERVATIONS 

The observations shown by bars in Fig. 1 were made 
by the ISEE-1 and -2 and Voyager- 1 and -2 plasma wave 
instruments. They show the ranges of rms Langmuir 
fields measured in the instrumental frequency channel 
nearest u p in each case. It should be noted that mea- 
sured Langmuir fields tend to be underestimates for 
three reasons: First, the waves are narrowband, so if 
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Figure 1. Electric fields of Langmuir waves observed 
in the foreshocks of Venus, Earth, Mars, Jupiter, Sat- 
urn, Uranus, and Neptune are shown vs. heliospheric 
distance. The theoretical upper bound (5) is shown as 
a sloping line for £ = 0.35. ISEE and Voyager obser- 
vations are shown as solid bars, while upper bounds 
obtained by other spacecraft are shown as dots. 
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Galileo and the Pioneer Venus Orbiter [Hospodarsky et 
al ., 1994, and references therein], much lower than (5)’s 
prediction of ~ 40 mV m -1 for 7 ^# = 10 _ 3 w p . Ulysses 
fast envelope sampler observations have identified large 
spikes (lasting only 1-2 ms or 1 - 2 samples) super- 
posed on broader Langmuir wave packets in the Jovian 
foreshock [Thiessen and Kellogg, 1993]. The spikes are 
not related to the electrostatic decay process (see be- 
low). However, the broader wave packets have maxi- 
mum levels of ~ 5 mV m -1 , consistent with the bound 
shown in Fig. 1. IMP 6 observations yielded a maximum 
field of 70 mV m -1 in the Earth’s foreshock [Filbert and 
Kellogg, 1979], above the predicted bound (for nominal 
parameters) but below 3£xo • 

4. DISCUSSION 

We have shown that the theoretical electrostatic de- 
cay threshold Elq is consistent with highest observed 
Langmuir fields in planetary foreshocks. 

Several assumptions were made in Sec. 2 to obtain 
the scaling of Elq with heliospheric distance r, yielding 
a threshold energy density proportional to the plasma’s 
thermal energy density. It could be argued that any 
process which yields a maximum field proportional to 
NT e would also be a candidate to explain the observed 
scaling with r. However, any alternative must also ac- 
count for the quantitative values of the maximum fields. 

One alternative mechanism is wave collapse, which 
yield maximum attainable fields (at the end of collapse) 

E m » ( 2Nk B TJeo) 1/ 2 , ( 6 ) 

(e.g., Robinson and Newman [1989]). However, ( 6 ) 
yields an upper bound ~ 60 times higher than (5) for 
the same parameters and 7 l 1 = 10 - 3 u; p — far above 
anything observed. Nonetheless, fields of this order are 
conceivable, provided they occur so rarely as to have 
escaped detection by the Voyager and Ulysses instru- 
ments. Robinson and Newman [1989] presented evi- 
dence that coherent wave packets do collapse occasion- 
ally in the Earth’s electron foreshock, starting from a 
much lower background level of random-phase waves 
(with the collapse threshold near the rms wave level) . It 
is consistent that the upper bound to the latter waves is 
set by the electrostatic decay threshold, with occasional 
wave collapse occurring. In this event, very few observa- 
tions would be expected to detect extremely high fields 
between the collapse threshold and the boundary ( 6 ), 
because collapse is very fast in its latter stages when 
the fields are highest. Hence, the electrostatic decay 
threshold would still represent an effective cutoff ex- 
cept in very extensive sets of statistics. It is worth not- 
ing in this context that Cairns and Robinson [1992a, b] 
were unable to find any evidence for wave collapse in 
the Jovian foreshock; instead, they found strong evi- 
dence of electrostatic decay, as did Thiessen and Kellogg 
[1993]. This conclusion has since been complemented 


by similar results from independent studies of type III 
source regions [Robinson et al., 1993a; Gurnett et al., 
1993; Cairns and Robinson, 1994]. Thiessen and Kel- 
logg [1993] also showed that most of the spike events 
observed by Ulysses in the Jovian foreshock are quan- 
titatively inconsistent with strong turbulence theory. 

As in foreshocks, Langmuir waves in type III radio 
sources are generated by beams and it is worth asking 
whether the above scalings also apply to them. Signifi- 
cantly, earlier studies of type III sources implied a varia- 
tion of the peak fields as r -1-4±0-5 for r < 1 AU [Gurnett 
et al., 1980], consistent with the present work. Indeed, 
observations where r < 1 AU imply that NT e falls more 
rapidly than for r > 1 AU, with NT e ~ r*" 1-3 * 0-1 , so the 
agreement is even better than for £ = 0.35. Using (2), 
we find Elq = (2 — 8 ) mV m -1 at 1 AU for typical type 
III parameters T e = (1 — 2) x 10 5 K, — (10 — 30) V, 
u>p — (1 — 2) x 10 5 s -1 , A vb/vb = 0.3, 7 = 1.3 — 2.0, and 
7 L < = 10 s -1 [Robinson et al., 1993a, b]. The highest of 
these upper bounds agrees well with the maximum ob- 
served value of 7 mV m -1 [Gurnett et al., 1980]. Satu- 
ration via electrostatic decay can thus explain the upper 
bounds of both foreshock and type-III Langmuir fields. 

Theoretically, the exponent £ in (4) for the radial 
variation of T e can be related to the electron ratio of 
specific heats 7 e by £ = 2(7 e — 1) [Scudder and Olbert, 
1979; Sittler and Scudder, 1980]. The latter authors 
fitted Mariner and Voyager-2 data from 0.45 to 4.76 AU 
and found very good agreement with (4) for £ = 0.35 ± 
0.06, although Phillips et al. [1993] obtained £ = 0.70 
and 0.24 by analysing Ulysses data from 1.2 to 3.8 AU 
and 3.8 to 5.3 AU, respectively. Alternatively, fitting 
(5) to the maximum wave fields in Fig. 1, not including 
Venus, leads to the independent estimate £ = 0.6 ± 0.5. 
This range is consistent with the electron observations 
but, except for ruling out adiabatic cooling with £ = 
4/3, does not discriminate between the other theoretical 
and experimental values of £ proposed previously. 

In conclusion, from the above analysis and discussion 
of foreshock and type III Langmuir wave fields, and 
from the previous work mentioned in both contexts, the 
electrostatic decay is chiefly responsible for saturating 
the highest levels of beam-driven Langmuir growth in 
both cases. 
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